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Abstract. We characterize correlations in polarization and time of photon pairs 
generated from a cold cloud of ®^Rb atoms via a four-wave mixing process in a cascade 
level scheme. Quantum state tomography reveals entangled polarization states of 
high purity for each of the decay paths through two different intermediate hyperfine 
levels. When allowing both decay paths, we observe quantum beats in time-resolved 
correlation measurements. 
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1. Introduction 

Time-correlated and entangled photon pairs have been an important resource for a wide 
range of quantum optics experiments, ranging from fundamental tests |1] to applications 
in quantum communication, cryptography, teleportation and computation j2]. 

The hrst sources of correlated photon pairs were based on a cascade decay in neutral 
atoms mil]. The cascade imposes a time correlation, and with an appropriate choice 
of the geometry and intermediate states it is possible to observe a strong non-classical 
correlation in the polarization of the photons. These sources are rarely used today in 
quantum optics experiments for their limited generation rates, and have been superseded 
by schemes based on three-wave mixing in non-linear optical crystals, or four-wave 
mixing (FWM) in optical hbers [5]. To obtain a narrow optical bandwidth, it is possible 
to use near-resonant transitions in atoms to provide the large third order nonlinear 
susceptibility for efficient FWM, leading to photon pairs with a central wavelength 
matching those of the transitions involved [HI [71IH] • These sources are a hybrid between 
the traditional atomic cascade approaches, and those based on three- or four-wave 
mixing in solids in the sense that they deliver both a useful pair rate collected into 
single mode optical hbers, and exhibit interesting temporal correlations. 

In this paper, we present the characterization of a source of time-correlated and 
polarization-entangled photon pairs based on four-wave mixing in a cold cloud of 
®^Rb atoms. The involved atomic levels, selected by the choice of pump and target 
wavelengths, form a cascade decay scheme, providing an asymmetrical time correlation 
similar to the one from the cascade decay of single atoms. 

2. Polarization entanglement 

Entanglement between photons can be established in several degrees of freedom like time 
bins, polarization, and orbital angular momentum [9], with polarization entanglement 
having been extensively studied due to its robustness and the availability of very stable 
optical elements for manipulating and detecting the polarization of single photons |T0] . 

Observation of photon polarization entanglement begin with early photon pair 
sources based on cascade decays in atomic beams m, followed by spontaneous 
parametric down conversion in nonlinear optical crystals na, cold |H] and hot PI 
atomic vapors, and recently also in cascade emission from quantum dots jll]. In this 
paper, we characterize the correlation polarization properties of nearly Fourier-limited 
photon pairs generated from the cold cloud of ®^Rb [IHl [16]. 

3. Experimental setup 

The experimental setup is similar to our earlier work but uses a collinear beam 

geometry (see £gure[I]). An ensemble of cold ®^Rb atoms is prepared with a magneto¬ 
optical trap (MOT) of optical density OD ^ 32 for light resonant to the 5Si/2, F = 2 
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Figure 1. Level scheme for four-wave mixing in ®^Rb, and schematic of the 
experiment. A first interference filter (IFi) combines the two pump beams in a co¬ 
propagating geometry inside the cloud, a second one (IF 2 ) separates the signal and 
idler photons from residual pump light. The pump polarizations can be freely chosen 
with polarizers (P) and quarter wave plates (q). A stack of quarter wave plate, half 
wave plate (h), and polarizing beam splitter (PBS) in each collection mode can select 
any arbitrary polarization. A solid etalon (E) can be used to select light from only one 
for the decay paths X and Y. Di, Ds: Avalanche Photodetectors. 


—^ 5 P 3 / 2 , F = 3 transition. The atoms are excited from 5Si/2,F = 2 to 5 D 3 / 2 , F = 3 
via a two-photon transition, with a two-photon detnning of ~ 5 MHz. Pump beams of 
wavelength 795 nm and 762 nm overlap in a co-propagating geometry inside the cloud. 
The 795 nm pump is red detuned by 30 MHz from the intermediate level 5Pi/2, T" = 2 to 
reduce the incoherent scattering rate. From the 5 D 3 / 2 , F = 3 excited level, atoms can 
decay through several paths. We select “signal” photons around 776 nm, and “idler” 
photons around 780 nm with interference hlters of 3nm FWHM bandwidth. Within 
this bandwidth, two decays can be observed (figure H] top right): decay X through the 
hyperfine level 5 P 3 / 2 , F = 3, and Y going through 5 P 3 / 2 , F = 2. 

Energy conservation and phase matching results in the generation of signal and 
idler photon pairs from both decay paths with a frequency difference of 5=266 MHz 
corresponding to the hyperhne splitting of the intermediate level. The generated photons 
are collected into single-mode fibers and detected by avalanche photodetectors (quantum 
efficiency ~40%, jitter time ~ 1 ns). In the experiment, we cycle between a 150/is long 
cooling period with the MOT turned on, and a 10 /rs long period for pair generation. 
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Figure 2. Tomographic reconstruction of the polarization state p (real part) for 
biphotons generated via decay path X (on the left), and Y (right), for pump modes 
set to orthogonal circular polarizations. The imaginary part of all matrix elements is 
below 0.09 and is not shown. 


4. Polarization state tomography 

We investigate the polarization state of the photon pairs for decay paths X and 
Y independently. We select light in the signal mode (see hgnre [T]) nsing a 2 cm 
long solid fnsed silica etalon with a transmission bandwidth of 52 MHz (FWHM). 
The etalon is temperatnre-tnned to match to the resonance freqnency of either the 
5 P 3 / 2 , F = 3 —>■ 5 D 3 / 2 , F = 3 or 5 P 3 / 2 , F = 2 — 5 D 3 / 2 , F = 3 transition; its 
temperatnre is stabilized to within 1 mK to minimize freqnency drifts. 

We completely characterize the polarization state of photon pairs via qnantnm 
state tomography ra by projective detections of individnal photons in a combination 
of linear and circular polarizations. For this, we insert quarter and half wave plates (q, 
h) followed by beam splitter cubes (PBS) in the signal and idler modes. Figure [2] shows 
the real part of the reconstructed biphoton states px and py corresponding to decay 
paths X and Y. The imaginary parts of all elements are smaller than ±0.09h The 
strong off-diagonal elements {LR, RL) signify the polarization entanglement. From the 
reconstructed matrices, we can extract typical entanglement measures for the state; we 
evaluate the concurrence C [18], and the entanglement of formation E [19]. Furthermore, 
we can also determine the purity P = Pr[p\ y] of the biphoton state for each decay path. 
The values of these indicators are given in table [H 

Despite the fact that the atomic ensemble is not prepared in a particular Zeeman 
sublevel, the polarization states for photon pairs from both decay paths show a 
remarkably high purity. This is compatible with theoretical models presented in [201 |2T] , 
which we briefly summarize here. 

In a cascade decay, polarization entanglement arises from indistinguishable decay 
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X 

Y 

Purity P 

0.92±0.02 

0.96±0.03 

Concurrence C 

0.89±0.01 

0.94±0.01 

Entanglement of formation E 

0.85±0.03 

0.98±0.01 


Table 1. Entanglement indicators for reconstructed states px,Y- The uncertainties 
reflect propagated Poissonian counting statistics of contributing coincidence events. 


paths, in our case provided by sufficiently degenerate Zeeman states of each hyperhne 
level. With the quantization axis along the beam propagation direction of all modes, 
we only drive transitions with Amp = ±1 with orthogonally circularly polarized 
pump beams. In parametric processes [22], the quantum state of the medium remains 
unchanged through the interaction [23]. Further, rotational symmetry of the atomic 
cloud in beam propagation direction implies angular momentum conservation. Along 
with the angular momentum selection rules, this limits the possible polarizations of the 
generated signal-idler photon pairs to \LR) and \RL). Since the process is coherent and 
\LR) and \RL) are indistinguishable otherwise, the resulting state IV’) of a target mode 
photon pair can be written as 

\'ijj) = aQ\LR) + exp(i0o) ai\RL). (1) 


The probability amplitudes ao,i and the phase 0o can be derived using a model based 
on the relative transition strength between different Zeeman sublevels [201121] as 


®o,i 


X. 


as,ai 


^ {Xag,aj) 
Q5,Q:/=ibl 


( 2 ) 


where asg are the helicities of the signal and idler photons, and Xas,ai is the product 
of relevant Clebsh-Gordan coefficients m that couple the individual {mp) states of the 
different hyperhne levels involved in the four-wave mixing process, and 

P 

™ _ \ ^ ^F^,l,Fe ^Fg,l,Fj^ 

Xas,aj / ^ — — l '~'mp — l,l,mp '^mp—as,as,‘mp '^mp,—ai,mp—aj ) 

mp=-Fg 

where Fg,b,e,d = 2, 2, 3, 3 corresponding to the respective total angular momentum F of 
the participating levels. From ([2]), we obtain the expected state I'lpx) ~ 0.55\LR) — 
0.83\RL) for the decay path X. The reconstructed state px matches the expected one 
with a hdelity of 94±1%. For the decay path F, the model predicts state iV^y) ~ 
0.92|Li?) — 0.39|i?L), which agrees with py with a hdelity of 93±1%. 


5. Transition strength of different decay paths 

Both decay paths exhibit diherent decay time constants due to diherent transition 
strengths for the decays. The transition 5 P 3 / 2 , F = 3 to 5Si/2, F = 2 is 2.8 times 
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Figure 3. Coincidence events as a function of the detection time difference At 
between signal and idler photon detection for the decay path X and F, selected with 
a temperature tuned etalon. Integration times were 7 and 14 minutes, respectively. 


stronger than from 5 P 3 / 2 , F = 2 [25]. This results in a higher optical density (OD) 
for the F = 3 transition [26], and should lead to a faster decay via path X m 
To experimentally investigate this, we perform separate time correlation measurements 
between the detection of signal and idler photons for each decay path. The histogram 
of coincidence events as a function of time delay At between the detection of signal and 
idler photons sampled into 1 ns wide time bins is shown in hgure [3j The solid line in 
both cases shows a £t to a heuristic model, inspired by the joint effect of a decay time 
corresponding to the two-photon transition [T5[ [28] , and the effect of a finite ring-down 
time of the filter etalon: 


G 


( 2 ) 

x,v 


(At) 


Go exp(At/Tr) for At < 0 

Go exp(—Af/rx,y) for At > 0 . 


(4) 


For path X, we obtain a decay constant rx=5.6±0.1ns for an idler photon heralded 
by a signal photon. In the same way, with the etalon tuned to transmit the resonance 
frequency of the 5 P 3 / 2 , F = 2—)- 5 D 3 / 2 , F = 3 transition for path Y, we find ry= 
13.1±0.2ns. Both decay constants tx,y are shorter than the spontaneous decay time 
Tgp = 27 ns of the 5 P 3/2 level of a single Rubidium atom in free space due to the 
collective enhancement effects observed in an optically thick atomic ensemble . 

The rise time = 3.1 ± 0.3 ns for decay path X and = 3.3 ± 0.4 ns for R is a 
consequence of the finite response time of the etalon. Both values are compatible with 
the value of 3.0 ± 0.1ns obtained in an independent characterization of the etalon. 


6. Quantum beats 

Without the etalon (see figured]), the decay paths X and Y cannot be distinguished by 
wideband photodetectors. Consequently, the energy difference between photons from 
the two paths leads to a modulation of the time correlation function between signal 
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At (ns) 

Figure 4. Coincidences as a function of the time delay between the detection of 
signal and idler photons, with no etalon in the signal mode. Pump 1 and 2 are 
set to orthogonal linear polarization H and V, respectively, and signal and idler 
are projected onto the polarization states \L) and (0.7 + 0.57z)|i7) + 0.41 j \V). The 
observed modulation (“quantum beat”) is associated with the hyperfine splitting of 
266 MHz between F = 3 and F = 2 of the 5 P 3/2 level. To resolve the oscillations with 
high contrast, avalanche photodetectors with a low time jitter (sidOps) were used for 
this measurement. Due to the lower quantum efficiency of these detectors, the total 
acquisition time is 5 hours. 


and idler photodetection, as shown in Figure 01 This, and other similar phenomena, is 
known as quantum beats: it was predicted at an early stage of quantum physics |29], and 
hrst experimentally observed in pulsed optical excitation of atoms with two upper states 
decaying to the same ground state [30l [3T]. Quantum beats have also been observed in 
cascade decay systems of dilute atomic vapours [32], dense thermal atomic vapours [33] . 
and for single ions |34j . 

In our case, the beat frequency 5/27r=266 MHz is equal to the energy difference 
between the hyperhne levels 5 P 3 / 2 , F = 3 and F = 2. The measurements shown in 
Figured] are performed with polarization of pump 1 and 2 set to H and V, respectively, 
and polarizations of signal and idler modes are set to observe a large interference 
contrast. 

To model the interference between the different decay paths, we take into account 
the different coherence time of the emitted photon pairs, characterized by the time 
constants tx and Ty, and the relative difference in the generation amplitude R, and 
phase 0. These last two terms can be calculated via (E]). We express the probability 
amplitudes cx,y for paths X and V as function of the detection time difference At 
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Figure 5. Coincidence rate as a function of time delay between the detection of signal 
and idler photons for different choices of polarization of signal and idler photons, (a) 
The beats are damped by choosing the appropriate polarizations due to suppression 
of coincidences from decay path Y. R — 2.86 * 10“^ and f = tt, total acquisition 
time 9 minutes, (b) and (c): different polarization projections change the phase of 
the oscillation. In this two cases, the relative phase difference is tt. (b) R = 1.43 and 
(/) = 0, total acquisition time 35 minutes; (c) R — 0.5 and (j) = n, total acquisition time 
28 minutes. 


between signal and idler photons, 

cx(At) = and cy(At) = , (5) 


which interfere to a joint detection probability 

= |cx + Cy|^ 


= 0(At)G^ 


At _At At(Ty +Ty) 

e + R^e w + 2Re cos {6 At + 0) 


( 6 ) 


Using the measured coherence times tx and ry and the values of R and 0 calculated 
from the interaction strengths of the transitions, we £t the experimental data in hgure H] 
using (E]) with only Go and an accidental count rate as free parameters, and find a good 
agreement with this relatively simple model. 

Different interaction strengths for the different polarizations in the participating 
levels allow control of the relative amplitude and phase of the possible decay paths. 
We can observe the dependence of the amplitude of the oscillation on the polarization 
settings and compare it with the expected values calculated from the interaction 
strengths of the transitions. In figure [5] we present three different cases, all fitted in 
a similar way as for hgure ID Of particular interest is the case where the beats are 
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almost entirely suppressed [figure [Sl[a)], an indication that most of the photon pairs 
observed are generated by the X decay. Figure [5](b) and (c) show the situation for 
polarization selections that lead to quantum beats with a high contrast, but opposite 
phases. 

This demonstrates that it is possible to select one frequency component only by an 
appropriate choice of polarizations, without using an etalon. However, it is not possible 
to select only photon pairs from the Y decay in a similar manner due to the relative 
weakness of the transitions involving the 5 P 3 / 2 , F = 2 level. 


7. Conclusion 

In summary, we have characterized the polarization entangled state of photon pairs 
from a cold cloud of atoms by performing quantum state tomography, individually 
for two decay paths of the cascade. We hnd that the resulting polarization-entangled 
states for both decay paths are not maximally entangled, but reasonably close to it. 
This is compatible with a model combining the transition strengths between different 
participating intermediate states in the four-wave mixing process. We observe high- 
contrast quantum beats in a time correlation measurement between the generated 
photon pairs. The contrast and the initial phase of beats can be controlled with the 
choice of polarization of pumps and projective measurements on the generated photons. 
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